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We study equilibrium shapes, stability and possible bifurcation diagrams of fluids in higher dimen- 
sions, held together by either surface tension or self-gravity. We consider the equilibrium shape and 
stability problem of self-gravitating spheroids, establishing the formalism to generalize the MacLau- 
rin sequence to higher dimensions. We show that such simple models, of interest on their own, also 
provide accurate descriptions of their general relativistic relatives with event horizons. The examples 
worked out here hint at some model-independent dynamics, and thus at some universality: smooth 
objects seem always to be well described by both "replicas" (either self-gravity or surface tension). 
As an example, we exhibit an instability afflicting self-gravitating (Newtonian) fluid cylinders. This 
instability is the exact analogue, within Newtonian gravity, of the Gregory-Laflamme instability 
in general relativity. Another example considered is a self-gravitating Newtonian torus made of a 
homogeneous incompressible fluid. We recover the features of the black ring in general relativity. 
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INTRODUCTION 



The shape and stability of liquid drops has been the subject of intense study for more than 200 years, with exciting 
revivals occurring more or less periodically. In the early days of the topic, the hope was that centimeter-sized drops 
held together by surface tension could be good models for giant liquid masses held together by self-gravitation. 
Liquid drops, which could be controlled in the laboratory, should model self-gravitating stars where gravity would 
replace surface tension. A great deal of what is now known about (rotating) liquid drops is due to the work of J. 
\Q ' Plateau |lj who, despite suffering from blindness throughout a good part of his scientific life, conducted many rigorous 
experimental investigations on the subject (we refer the reader to pj for details). Starting with a non-rotating drop of 
spherical shape, Plateau and collaborators progressively increased the rate of rotation (angular velocity, assuming the 
drop rotates approximately as a rigid body). Plateau found that as the rotation rate increased, the drop progressed 
through a sequence of shapes which evolved from axisymmetric for slow rotation, to ellipsoidal and two-lobed and 
finally toroidal at very large rotation. These experimental works have been repeated with basically the same outcome, 
. even in gravity- free environments (see also Q)- Theoretical analysis of Plateau's experiments were started by 
Poincare Q, who proved that there are two very different possible solutions to the equilibrium problem: one is 
• i-h , axisymmetric like disks or torii and the other is comprised of asymmetric shapes with two lobes. The analysis was 
completed by Chandrasekhar improved upon by Brown and Scriven |(J (see also 0,0) and is summarized in Figure 
^ taken from 6]. The Figure shows the (dimensionless, as defined by Brown and Scriven, see caption in Fig. ^) 
angular velocity as a function of the angular momentum of the drop, and it displays the stability and bifurcation points 
of the evolution. Let's focus first on the evolution of axisymmetric figures. Starting from zero angular momentum, 
where the shape is a spherical surface, one increases the angular momentum. The resulting equilibrium shape evolves 
from a sphere through oblate shapes to biconcave shapes (see also |5j). There are no (axisymmetric) equilibrium 
figures with SI > Sl max ~ 0.753. The axisymmetric family does not terminate at SI = Sl max , instead it bends back to 
lower angular velocity, and there are now two axisymmetric figures of equilibrium. The steepest descending branches 
are torii 0, |(| ■ 

Stability calculations show that all shapes of the axisymmetric family are stable against axisymmetric perturba- 
tions, however a real drop is also perturbed by non-axisymmetric modes. The stability calculations show that the drop 
becomes unstable at SI = S7// ~ 0.5599 [||, where the axisymmetric family is neutrally stable to a two-lobed perturba- 
tion. This point marks the bifurcation of a two-lobed family (therefore a non-axisymmetric family) . All axisymmetric 
shapes with SI > SI// are unstable against two-lobed perturbations, but there are additional points of neutral stability, 
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FIG. 1: The bifurcation and stability diagram for rotating liquid drops with surface tension (self-gravity is neglected). Shapes 
are stable along solid branches, unstable along broken ones. The plot depicts dimensionless angular velocity as function 

of dimensionless angular momentum. The dimensionless angular velocity is Q, — \J where all hatted quantities are 

dimensionfull, R is the radius of a sphere with the same volume as that of the drop, p and T are the density and surface tension 

of the liquid. Likewise, the dimensionless angular momentum J is defined as J = J/4TR 2 \JiTR 2 pJ. Although not shown in 

the plot, the axisymmetric curve bends back to lower angular velocities after it has passed the four-lobed family neutral point. 
Note also that the maximum allowed angular momentum for stable shapes is denoted by Ju in the diagram. Taken from |g. 



which are also bifurcation points: the three- and four-lobed perturbations are shown in the diagram Q but there are 
others |8(. The first, most important non-axisymmetric surface, the two-lobed family, branches off at fi = and 
is initially (near the bifurcation point) stable. In Figure [21 we show some typical members of that sequence (see [lOj, 
where this figure was taken from, for further details). It turns out that the evolution of self-gravitating fluids, i.e., 
stars, follows a very similar diagram, as shown in Figure Finding the equilibrium configurations of self-gravitating 
rotating homogeneous fluids is an extremely important problem, which has met intense activity ever since the first 
discussions by Newton and others (see [TJ and also the excellent and concise introduction in [12), which was revis- 
ited some 50 years ago by Chandrasekhar and co-workers |l2]. In three spatial dimensions, the most basic figures 
of equilibrium of rotating bodies are the MacLaurin and Jacobian sequences j7fij . The MacLaurin sequence is a 
sequence of oblate spheroids along which the eccentricity of meridional sections increases from zero to one. Jacobi 
showed that for each value of f2 2 , the square of the angular velocity of rotation, there are actually two permissible 
figures of equilibrium. The other is a tri-axial ellipsoid, belonging to the Jacobi sequence, which branches off from 
the MacLaurin sequence. This is the first of several points of bifurcation which distinguish the permissible sequences 
of figures of equilibrium of uniformly rotating homogeneous masses. If one starts from a non-rotating homogeneous 
fluid ellipsoid in the MacLaurin sequence and increases the angular momentum of the body, one will first come across 
the Jacobi branch, where the ellipsoid becomes unstable with regard to the first non-axisymmetric perturbation [l^ 
(which are the self-gravitating analogues of the shapes in Figure Continuing, still along the MacLaurin sequence 
and increasing angular momentum, one next comes to the bifurcation point of other non-axisymmetric figures: the 
"triangle", the "square" and "ammonite" sequences [l2l[T^ |. One finally encounters the first axisymmetric sequence. 
The bodies of this se quence p inch together at the centre and eventually form Dyson's anchor-rings, or torus- like con- 
figuration |15L Ho . Ii7l flSl. fl9l | . It is apparent that the bifurcation diagram for self-gravitating fluids, Figure is very 
similar to the one for non-gravitating fluids with surface tension, shown in Figure^ Thus, whenever the treatment of 
self-gravitation proves too difficult, the surface tension model is a good guide. A realistic star is neither homogeneous 
nor incompressible and it certainly isn't held together by surface tension. However, it seems logical to build our under- 
standing using simpler models. We have gone a long way with simpler related models: another well-known example, 
where the fluid drop model proved extremely useful, is Bohr and Wheeler's po| proposal to describe nuclear fission as 
the rupture of a rotating charged liquid drop, where now the surface tension plays the role of nuclear forces. A review 
with a detailed and unified discussion of charged nuclei, astronomical bodies and fluid drops can be found in (a 




FIG. 2: This plot shows a sequence of members of the two-lobed family of drop shapes for several values of angular velocity. The 
two-lobed drops seem to approach a limiting surface consisting of two spheres touching each other in one point as Q approaches 
0. The first shape, with Q = 0.5599, corresponds to the first (two-lobed) bifurcation point in Fig. Q For SI = 0.2885 for 
instance, the figure is peanut shaped, and we can situate the corresponding family in Fig. Q Taken from (l0| . 
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FIG. 3: The bifurcation diagram for rotating self-gravitating, incompressible homogeneous stars. Here we show the angular 
velocity as function of the angular momentum for several possible families. The data and notation are those of Eriguchi and 
Hachisu 0. The line A — A 1 is the MacLaurin sequence, B — B' the Jacobi sequence, C — C triangle sequence, D — D' 
the square sequence, E — E' ammonite sequence, F — F' the one-ring sequence, G — G' the two-ring sequence, H — H' the 
pear-shaped sequence, I — I' the dumb-bell sequence, J — J' the Darwin sequence and finally K — K' is the two point-mass 
sequence. Note the similarities with Fig. See [l4| for further details. 



concise and complete account of the history of this subject, along with some good quality figures can be also found in 
|10|). Liquid drops can model self-gravitating systems. Can they model extremely compact, highly relativistic objects 
like black holes? The answer, quite surprisingly, seems to be yes |2^. Regarding the event horizon as a kind of fluid 
membrane is a position adopted in the past [23. A simple way to see that this is a natural viewpoint, is to take the 
first law of black hole mechanics |24j which describes how a black hole (we will take for simplicity uncharged, static 
objects), characterized by a mass M, horizon area A, and temperature T = 1/(32ttM), evolves when we throw an 
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infinitesimal amount of matter into it: 

dM = TdA . (1) 

We now know that the correct interpretation of this law is that black holes radiate, and therefore can be assigned a 
temperature [2f|,[2(|. The second law of black hole mechanics dA > is then just the second law of thermodynamics. 
But we can also argue that Equation (JJJ can be looked at as a law for fluids, with T being an effective surface 
tension 1271 . The first works in black hole mechanics actually considered T as a surface tension (see the work by 
Smarr |28| and references therein), which is rather intuitive: the potential energy for fluids, associated with the 
storage of energy at the surface, is indeed proportional to the area. Later, Thome and co-workers [23J developed 
the "membrane paradigm" for black holes, envisioning the event horizon as a kind of membrane with well-defined 
mechanical, electrical and magnetic properties. Not only is this a simple picture of a black hole, it is also useful for 
calculations and understanding what black holes are really like. This analogy was put on a firmer ground by Parikh 
and Wilczek p9l |. who provided an action formulation of this membrane picture. There are other instances where 
a membrane behavior seems evident: Eardley and Giddings [3(j, studying high-energy black hole collisions found a 
soap bubble- like law for the process, while many modern interpretations of black hole entropy and gravity "freeze" 
the degrees of freedom in a lower dimensional space, in what is known as holography |3l| . In |32| it was shown that a 
"membrane" approach works surprisingly well, yielding precisely the same results as the AdS/CFT correspondence. 
There have also been attempts to work the other way around: computing liquid surface tension from the (analog) 
black hole entropy [3^. As a last example of a successful use of the membrane picture, Cardoso and Dias have 
recently been able to mimic many of the aspects of a gravitational instability by using a simple fluid analogy. We will 
comment on this work later on. 

Having said this, it strikes one as very odd that the study of simple fluid models has not been extended to higher 
dimensions. Extra dimensions, besides the four usual ones, seem to be needed for consistency in all or many of 
candidate theories of unification. We have in mind not only string theory |34j where the extra dimensions are 
compactified to a very small scale, but also recent theories with large |35l | or warped extra dimensions |36j |. Therefore, 
an understanding of higher dimensional gravitating objects is imperative. However, even though there are several 
known exact solutions to Einstein equations in more than four dimensions describing black objects (with event 
horizons), nothing or very little is known about the simpler Newtonian fluid equilibrium configurations, with or 
without self-gravity. 

For instance, higher dimensional black hole solutions with spherical topology have long been known |37l l38| , which 
generalize the four-dimensional Kcrr-solution to higher dimensions. In four spacctimc dimensions one can show that 
the topology of a black hole must be a sphere, but that's no longer the case in higher dimensions. A simple example 
is to consider extra flat dimensions, for instance Schwarzschild x R, which is still a solution to the vacuum Einstein 
equations. These extended solutions can then have different topologies, like a cylinder, and are called back strings 
|39|. Yet another example is the five-dimensional black ring |4(J,|4l|], having the topology of a torus. 

It is curious to note that we are in possession of several exact solutions to the full Einstein equations and yet 
nothing is known about the simpler case of self-gravitating Newtonian fluids or even fluids without gravity but with 
surface tension. To be specific, we found a solution as "strange" as a black ring in five dimensions, yet we do not know 
how a rotating, self-gravitating Newtonian fluid behaves as a function of the rotation rate, for higher dimensional 
spacetimes. These simpler objects should be much easier to study, while carrying precious information about their 
more compact cousins with event horizons. There are several topics that could and should be addressed within this 
research line: 

(i) Do the fluid equations admit toroidal-like configurations, which may collapse to black rings for very high densities? 
Are these stable? 

(ii) Ultra-spinning black holes in spacetime dimensions larger than four are conjectured to be unstable [2^, [42] for 
large enough rotation rates. If the four-dimensional results for self-gravitating fluids or for fluids with surface tension 
carry over, this is to be expected. As we saw in the evolution diagrams, axisymmetric configurations are unstable for 
rotation rates larger than a given threshold, decaying to a two-lobe, asymmetric configuration. This suggests that 
the end-state of ultra-spinning black holes is a two lobe, asymmetric state like those in Figure |21 This possibility was 
already put forward in [42j, but the arguments presented here make this a much stronger possibility. These shapes 
would then presumably spin-down due to radiation emission. Can we use simple fluids with surface tension or with 
self-gravity to model this process? Can we estimate the rotation rate at which instability occurs? 

(iii) The black strings and p-branes are unstable against a gravitational sector of perturbations. Can we simulate 
this instability with fluids? A first step was taken in [22| , where it was shown that a simple fluid with surface tension 
can account for most of the properties of this instability. Can we add self-gravity to the model? Can we include 
rotational and charge effects? 

(iv) Can we build the evolution and bifurcation diagrams for fluids with surface tension (or with self-gravity) in 
higher dimensions? Can we deduce other solutions to the general relativistic equations using these simpler models? 
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Can we infer other properties of already known solutions? 

These are just some of the topics that could be addressed with simple Newtonian fluid models, but the list could 
easily be made much longer. The goal of the present work is to pave the way to some of these possible developments, 
but especially to show that these simple models do yield relevant information. We will attack the problem in its two 
versions: fluids held together by surface tension and by self-gravity. 

Plan of the paper 

This paper is divided in two parts, Part I in which we deal with fluids held by surface tension and Part II where 
we deal with fluids held together by self-gravity. 

We start with Section [n] where we review the meaning of angular momentum in higher dimensions and also the 
symmetries of a rotating body. We also specify the two most important cases (case I and II) we shall deal with 
throughout the paper. In Section IIIII we lay down the basic equations of hydrostatics and hydrodynamics in higher 
dimensions and in Section IIVI we write the formalism necessary to analyse small deviations for equilibrium. The 
formalism just laid down is then used in the next sections. In Section we study the equilibrium figure of rotating 
drops. Oscillation modes of non-rotating spherical drops are derived in Appendix lAl and constitutes a generalization 
to arbitrary number of dimensions of the well-known result in three spatial dimensions. In Section I VII we discuss 
equilibrium torus-like configurations of rotating fluids held together by surface tension, and show that these share 
common features with their vacuum general-relativistic counterpart, the black rings. In Section IVIII wc study the 
stability of the equilibrium figures discussed in Section ^ against small perturbations. We compute neutral modes 
of oscillation and the onset of instability which allows one to derive the bifurcation diagram for the case of one non- 
vanishing angular momentum. We finish t his p art with Section IVIIII by reviewing the Rayleigh-Plateau instability of 
a thin long cylinder in higher dimensions |22(, which could be a good model for the Gregory-Laflamme instability. 
We discuss possible effects of rotation and magnetic fields. 

In the second part we deal with self-gravitating fluids in higher dimensions. We start in Section llXI by considering 
an infinite cylinder of self-gravitating homogeneous fluid. We show that this cylinder is unstable against a classical 
instability, here termed Dyson-Chandrasekhar-Fcrmi (DCF) instability, which we demonstrate to be the exact analogue 
of the Gregory-Laflamme |3| instability of black strings. We discuss implications of this result for rotating and charged 
fluids and black strings. We then discuss in Section|X]toroidal equilibrium configurations of self-gravitating fluid (often 
termed Dyson rings) and we compare them to exact general relativistic toroidal- like configurations, the black ring 
solution. Finally, in Section lXll we write down the basic equations that generalize the MacLaurin and Jacobi sequence 
to higher dimensions, equations which could be used in future numerical work to reconstruct the whole bifurcation 
diagram. Oscillation modes of non-rotating globes are derived in Appendix [5] an d constitute a generalization to 
arbitrary number of dimensions of the well-known result by Kelvin in three spatial dimensions. We close in Section 
IXIII with a discussion of the results, implications and suggestions for future work. 

This work is by no means complete. Many important topics are left for future work, we discuss some of these in 
the concluding remarks. 

Notation 

There are many different symbols used throughout the text. For clarity we find it useful to make a short list of the 
most often used symbols. 

D — > total number of spatial dimensions. 

x l = (x a ,x A ) — > coordinates of a rotating body. We consider two particular cases: in case I there is a single 
rotation plane, and a = 1, 2; A = 3, D; in case II, D is odd and there are ^pL rotation planes. In this case 
o=l, ...D — 1; A = D. 

Cjj — > the area of the unit (D — l)-sphere, Cd = ^7-5-7- 

r (."2" ) 

m — > radial coordinate in cylindrical coordinates. 

lj — > Fourier transform variable. The time dependence of any field is written as e . In the case of a stable 
oscillation, it is related to the frequency by lu 2 = —a 2 . Note that pure imaginary values of uj mean a stable 
oscillation. 
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£1 — > angular velocity. 

B(x) — ► Gravitational potential. 

G — ► Newton's constant in I?— dimensions, defined by writing the gravitational potential as 

This is perhaps not the most common definition. A popular definition of Newton's constant is the one adopted 
by Myers and Perry [3^] and the relation between the two is derived in Appendix lEl 

T — > magnitude of surface tension, defined by P — Tn l i where P is the pressure in the interior of the body (we 
take the exterior pressure to be zero). We take n to be the outward normal to the surface the body. 

P — > Pressure. 

V — > Volume. 

p — » density. 

A — > wavelength of a given perturbation. 

fc — ► wavenumber = 2ir/\. 

P 3 = J v px l x 3 d D x — > Moment of inertia. 



Part I 

Fluids held by surface tension 



II. SYMMETRIES AND GEOMETRICAL DESCRIPTION OF A MULTI-DIMENSIONAL ROTATING 

BODY 



Let us consider a stationary rotating object in D-dimcnsional Euclidean space, with coordinates x l — [x , . 



')• 



The rotation group in D dimensions, SO(D), has D(D — l)/2 generators J 13 = Each generator J lJ corresponds 

to an 5*0(2) rotation, on the plane x % — x 3 . 

In a suitable coordinate frame J = (J lJ ) can be brought to the standard form (see for instance [38|) 



/ .h 
-Ji 



V 



J 2 
-J 2 o 



/ 



(3) 



This means that in this frame the rotation is a composition of a rotation in the plane x 1 — x 2 , a rotation in the plane 
x 3 — x 4 , and so on, without mixings. The rotation is then characterized by [D/2] generators (where [•] denotes the 
integer part), which are the generators of the Cartan subalgebra of SO(D) (namely, the maximal abelian subalgebra 
of SO(D)). In other words, in D dimensions there are 



N = 



(4) 



independent angular momenta. 

We will consider two particular cases: 

(I) Only one angular momentum is non- vanishing: J\ ^ 0, Ji>i = 0. 

(II) All angular momenta are non-vanishing, and they are equal: J\ 
D to be odd. 



Jj\[; in this case we also assume 



These will be referred to as case I and II respectively throughout the whole article. 
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A. Case I. One non-vanishing angular momentum 

In this case we are selecting a one-dimensional rotation subgroup 

SO(2) c SO(D) . (5) 

This means selecting a single rotation plane, x 1 —x 2 . We assume our object has the maximal symmetry compatible 
with this rotation: 

G = SO{2) x SO(D - 2) c SO(D) . (6) 
If we decompose the coordinates under G as 

x i = (x a ,x A ); a = 1, 2; A = 3, ...,D , (7) 

and wc define 



w = VVF+(^F, 

z = ^(x 3 ) 2 + • • • + {x D f , (8) 
the surface of our object depends only on m and z, and its equation can be written as 

z = f(w) . (9) 

i.e. 

H = *-/H = 0, (10) 

therefore the surface does not depend on the angle in the x 1 — x 2 plane, and it does not depend on the D — 3 angles 
among x 3 , . . . ,x D . We also define 

Notice that it must be <f> = — oo when z = for the surface to be smooth. 
Being w^ a — 2-, the normal to the surface is 

-0— ,— J (12) 

with modulus S^'S^'Sj — 1 + </> 2 , so the normal unit vector is 

n l = (-d> f- 1 • (13) 

Since the pressure is related to surface tension via P — Tn\, we will need a general relation to compute n\. Let us 
then compute n\ = n a a + n A A for this case I. We have 

<b' d> Id wd> 



a 1 D — 3 

Vi + 2 2 

B. Case II. All angular momenta coincident, and D odd. 

In this case the rotation group is broken to its Cartan subalgebra 

{SO{2))^ L c SO{D) . (15) 
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The rotation planes are 



~1 _ -.2 3 _ „4 D-2 _ D-l 



and the x 15 axis is left unchanged by the rotation[2j]. The decomposition of the coordinates is now 



x % = (x a ,x A ); a =1,...,D -1; A = L> . 



We define 



w = ^(x 1 ) 2 + {x 2 ) 2 + --- + {x D - 1 ) 2 , 

„ — „D 



The surface does not depend on the D — 2 angles among the x a . We also define 

Notice that it must be = — oo when z = for the surface to be smooth. 
Being zu M — — , the normal to the surface is 



£,<= {-</>— ,o) , 

with modulus ^S^Sj- = 1 + 2 , so the normal unit vector is 



n = -c 



,0 . 



Finally, 



1 d / w s - 2 c 



TO- - 2 (to + ^ 



(16) 

(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



III. HYDRODYNAMICAL AND HYDROSTATIC AL EQUATIONS 

Let X = (X 1 ) the position vector of a fluid element in the inertial frame, and x = (x l ) its position vector in the 
rotating frame, i.e. 



X = Rx 

where R = Ri or Rn according to the case under consideration and 



(23) 



/ cos fit sin fit 
— sin fit cos fit 



Ri = 



V 



\ 



/ cos fit sin fit 
— sin fit cos fit 



Rn = 



1/ 



cos Qt sin fit 
— sin fit cos fit 



V 



(24) 



1/ 



are D-dimensional rotation matrices in the cases I, II respectively. We introduce the velocity V and acceleration A by 

(25) 



X = RV 
X = RA 



where an overdot denotes time differentiation. 



9 



The hydrodynamical equation, in the inertial frame, is 



(with P pressure and p constant density of the fluid) which, applying gives 

A '—\h F - < 27 > 

Let us compute A. We have 

V = R T Ryi + x (28) 

and then 

A = R T RV + V = (R T R) 2 x + 2R T R± + x. (29) 

We have that the only non-vanishing components of R T R are 

(i? T i?) afc = -fte Qb (30) 

where e ab = e^ ab \ and e 12 = 1 in the case I, e 12 = e 34 = • • • = eC 23 " 2 )^- 1 ) = 1 in the case II. 
Therefore, in both cases ((R T R) 2 ) ab = tt 2 5 ab , and 

1 BP 

x a = ^— + fiV + 2tte ab x b (31) 

p ox a 

^ p dx A ^ ^ 

The second term in the right hand side of (|31|l is the centrifugal force and the third term is the generalized Coriolis 
force. At equilibrium x 1 — x l — 0, therefore 

§ - ' nV - 

9P 0. (33) 



dx A 



thus 



P = Po + ^pfl 2 w 2 . (34) 

This formula is valid both in case I and in case II, i.e., when there is one non- vanishing angular momentum and when 
all of them coincide (and D is odd). 



IV. THE VIRIAL THEOREM, SURFACE-ENERGY TENSOR AND SMALL PERTURBATIONS 

We would ultimately wish to assess stability of rotating drops, in order to generalize Fig. ^to higher dimensions. 
This calculation would also provide some hints as to whether or not one should expect to see instabilities in black 
objects in higher dimensions, as we stressed in the Introduction. To attack this problem we will follow Chandrasekhar's 
approach, to which we refer for further details. Alternative (numerical) approaches, perhaps more complete but 
more complex also, are described by Brown and Scriven Q and also Heine 0. 



A. The virial theorem in the inertial frame 



Consider a fluid of volume V with uniform density p, bounded by a surface S. The external pressure, for simplicity, 
is taken to be zero. Due to the existence of surface tension, the pressure P immediately adjacent to S, on the interior, 
is given in terms of the surface tension magnitude T by 



P = Tn\i . 



(35) 
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Multiplying (|26|l by and integrating over the volume V of the fluid we get 

±J Y PX*X^X = 2^-J Y X^X t (36) 

where X 1 ^ = | J v pX % X^d D X is the inertia tensor. Integrating by parts the second term on the R.H.S. of (|36|l we 
have 

X 3 ^-d D X = f X 3 PdS l -5 l] f Pd D X, (37) 



IV Jg Jy 

where dS l — n l dS and dS is the (D — l)-dimensional surface element on S. Using (|33|l we further have 

X i^-d D X = T [ X j n i i dS i -6 ij [ Pd D X. (38) 

V OX 1 Jg Jy 

Introducing this result in (|36|l one gets the virial theorem 



dt 

where 



^ / pX i X j d D X = 2T j + + 5 ij U, (39) 



n = / Pd D x, (40) 

Jv 

& ij = -T I X j n\dS l . (41) 



The quantity & 3 is the surface-energy tensor |5J- Conservation of angular momentum guarantees that & J is sym- 
metric. Indeed, 



d 
dt 



p(X l X ] - X ] X l )d D X = & 3 - & ]i = . (42) 

v 



B. The virial theorem in the rotating frame 

In a rotating frame of reference the equations of motion are written down in Il-i2[l . Multiplying those equations by 
x J , integrating over V and using the previous sequence of transformations we get 

^ / px a x j d D x = 2Z oj + e aj '+ S aj U + n 2 I aj + 2Sle ac f px c x 1 d D x, (43) 

dt Jy Jy 

4 I px A xH D x = 21 Aj + e Aj +S Aj Il, (44) 
at Jy 

where 

I a i = [ px a x j d D x. (45) 



v 



The tensor & J is symmetric in the rotating frame, too, as can be shown by transforming it with the rotation matrix 
R. 

We remind the reader that according to our conventions the indices i,j can take any values between 1 and D, but 
the indices a, b and A, B are restricted according to the situation we consider (case I or II), described in Equations 
(J7J) and i|17|) . In the following we always consider the rotating reference frame. 

C. Virial equations for small perturbations 



Having solved for the equilibrium equations and found the shape for a given rotating drop (this problem will be 
explicitly worked out in the next Section 0, we ask whether this is a stable equilibrium configuration or not. To 
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answer this we follow the usual procedure to assess stability: we perturb the equilibrium configuration. Suppose then 
that a drop in a state of equilibrium is slightly disturbed. All background quantities Qq are then changed by an 
amount 5Qq. The motions are described by a Lagrangian displacement of the form 

i(x)e^ , (46) 

where lj is to be determined from the problem. In the case of constant density, £ is divergence free. Define 



p^x j d D x , 



where the semi-colon indicates that a moment with respect to coordinate space is being evaluated. Define also 



Then the perturbation equations read p| 

uj 2 T a '' j - 2wfie a T cy = S6 aj + n 2 T ai + 5 aj SIl . 



2 r A 'i = 6G Aj + 8 Aj sn. 



The surface integrals 8& 3 have the form 
5& 3 -- 



T5 I x j n k k dS l = — T 



^n k k dS l + / x j 6(n k k )dS l + / x 3 n k k 5{dS l ) 



-T 



[ en k k dS l + [ x J 5{n k k )dS l 
J s Js 



s 

x 3 n 



k dx i 



dS l 



where we have used the general result (see for instance |5j) 



S(dS i ) = £ j .dS 4 



dx' 



dS 3 



and ■ = 0. Furthermore, we have (see 



n n —r \ n" -n J 

ox 1 I ox 1 



Equations (|4*9"|) . ijSTjl) can be decomposed as 



uj 2 T a ' b - 2Me ac T c < b 


= se ab ^ 


- n 2 r ab + 


tu 2 r A > B 


= 5& AB 


+ S AB SU 




= 5S aA - 


f- n 2 r aA , 


uj 2 T A ' a 


= 8G Aa . 





(47) 



(48) 



(49) 
(50) 



(51) 
(52) 

(53) 



(54) 



sn, (55) 

(56) 
(57) 
(58) 

These are the basic equations we will use later on to perform a stability analysis of rotating drops in higher dimensions. 
To complete the calculations we first need to compute the equilibrium figures of rotating drops and we this in the 
following. 



V. EQUILIBRIUM SHAPE OF A ROTATING DROP 
A. Case I. One non-vanishing angular momentum 
From equations {T3J , J3U , CIEl we find 

T 



1 



2„2 



P + ^ptt ™ 



1 d 



UJCi 



1 D - 3 



■ dm ^/l + cj? ^1 + 02 z 



(59) 
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thus 



\p q w 2 + \ P n 2 w A = -T-^L= + (D — 3)TJ + C (60) 



where C is a constant, and / is the integral 



/M=r — . (6i) 

7o fV)Ji + {&) 



If the drop encloses the origin (which we suppose for the time being), then va = must belong to the surface; being 
1(0) = 0, this means that C = and equation (|60fl reduces to 

^Po^ 2 + lpQ 2 ru 4 = -T + (D - 3)TI . (62) 



2 8' ^/TT? 

When D = 3 this equation can be solved analytically p|, but when D > 3 this is not possible, because equation 1621) 
cannot be expressed in the form ^f- = F(vj). 

We thus have to numerically integrate equation 1)59(1 . which we re- write as 

* + -IpOW = (-- — ^L= + -J=^) , (63) 

We can re-scale coordinates w — > aro and z — > az, where a will henceforth be the equatorial radius of the equilibrium 
figure. In this way all quantities are dimensionless and equation l|63|) is written as 

P a 2 Id wcj) 1 D-3 

— +4S a n7 — -— =+ , (64) 



with E a = gT . 

This equation can be integrated numerically, starting at w — where z = z(0) and z'(0) =0 = 0. In this sense we 
are solving an eigenvalue equation: given S Q , Pqcl/T and z(0) are determined by requiring that at the equator r = a 
we have z(a) = and z'(a) = — oo. The problem is highly non-linear, so the numerical procedure involved varying 
both P^ajT and z(0) simultaneously over a rectangular grid, in order to determine those values for which the outer 
boundary conditions (z(a) = and z'(a) — — oo) were satisfied. The accuracy gets worse for higher S and D but one 
usually has a precision of about 5% in our numerical results. Numerical results for D — 4 are shown in Figure 0] for 
some values of E s , which can be compared to the D = 3 results presented by Chandrasekhar in For S a = 0, the 
non-rotating case, we start by verifying numerically that the sphere is a solution, as it should, with P a/T = D — 1. 
The equilibrium figures are very similar and indeed the similarity seems to continue for larger D. The three different 
equilibrium shapes shown in Figure 0] all enclose the origin. As in the D = 3 case (see 0), we find that for E a larger 
than a certain critical value, the drop no longer enclosed the origin, giving rise to toroidal shapes. The condition to 
be satisfied at the critical E is that z = at r — 0. For D = 3 the critical value of E a is E a ~ 2.32911. For D > 3 
it is harder to determine the critical parameter, since there is no integral expression for the shape. Nevertheless we 
determine E cr jt ~ 2.3 for D = 4, 5 with a 5% uncertainty. An open important problem is to determine if a critical 
point exists for all D. If not, this might mean that the bifurcation curve^does not bend down, and therefore that the 
angular momentum is unbounded, which might be a good analogue for what happens with ultra-rotating |38j black 
holes. This may require more accurate numerical integrations, since we cannot trust our method for D > 5. 



B. Case II. All angular momenta coincident, and D odd. 



In case II we will find that the results are very similar to those for D = 3. From equations l|22|) . (|34|) . (|35|l we find 

1 Q 2 W 2 = T d_ ( W D ~^ 

2 P w w d - 2 dm I yrr? 



P + -ptfm 2 = --^zo^- ^=4= , (65) 



which is immediately integrated to give 
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FIG. 4: Typical (axi-symmetric) equilibrium figures of rotating drops for case I (one non- vanishing angular momentum), in 
D — 4. The (blue) dashed line corresponds to a spherical surface with E a = 0. The (magenta) dotted line corresponds to 
E a = 0.46 and the (black) solid line corresponds to the limit E a = 0.86. In all three cases the equatorial radius is unity. 



where C is a constant. 

Let's first suppose the drop encloses the origin. In this case, there are points where w = and therefore we must 
have C = 0. Equation (|(jfc>|) can then be written as 

At the equator, where tu = a we must have <fi — > — oo and then Ijbbfl implies that 

= 1 - E& , E & = . (68) 



(D-1)T 2(D + 1)T 

Equations Hb7[) - I|68|) are the same as equations (A4)-(A5) by Chandrasekhar So, figures of equilibrium enclosing 
the origin do not depend on the dimensionality of spacetime, if measured in terms of the dimensionless number E b 
(by this we mean of course that z(w, E b ) is independent of D; anyway, the relation between E b and Q depends on D, 
thus z(zu, fl) depends on D). By measuring w and 4> in units of a, we get 

t = ±= "(i-s» + ^ a ) (69) 

dm 



1 — vj 2 (1 — E5 + T,i,zu 

which can be written as 



z = /(tu) = / dw . (70) 

y/l- u; 2 (1 - E b + E feW 2 ) 2 



Transforming to a new variable 9 defined through 



w 1 = l-Atan 2 (0/2), A = V1 ± 2E " , (71) 

Efc 



the integral representation can be simplified to 



1 r ( „ 2AE b \ d9 

1 + ^- T -— -A ~ n = . ., = , (72) 



2E b% /A.7o V l + cos0/ ^! _ x 2 sin 2 
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where 



1 I 1 



Here, <j) is related to the value of w at the boundary by 



'1 -.4 tan 2 0/2 

The solution O can be expressed in terms of elliptic integrals 

dO 



as 



F(X, 



1 



y/l 



X 2 sin 2 i 



E(X, 



dOdl-x 2 sin^ 



2E,,VA 



(l-Z b A)F( x ^) + 2AX b E( X ,< 



sin 



X 2 sin 2 i 



1 + COS ( 



(73) 



(74) 



(75) 



(76) 



Some representative equilibrium solutions are displayed in Figure 03 for some values of E&. The shapes are very 
similar to their "case I" relatives. We show typical shapes for high values of E (these can be compared with the ones 
presented in Figure 1 of Chandrasekhar [5j). 

We have been dealing with drops that enclose the origin, and found that a closed expression for the shape is given 
by (|76|l . When does this sequence end? A necessary and sufficient condition for the drop to enclose the origin is that 
z = and £ = are simultaneously satisfied. Now, the condition £ = implies that <fi = </> m ax = 2 arctan A^ 1 / 2 . 
Substituting this value in l(7(3|) we get 



E max = 2.32911 . 

For values of E;, larger than the above, the figure no longer encloses the origin. 



(77) 



VI. EQUILIBRIUM TOROIDAL SHAPES 

In the general case of figures which do not enclose the origin (i.e.,zu — is not attained), we have to use equation 
(15911 and (|66|l for cases I and II respectively. Here we shall focus on case II, equation Ij66|l . since it is easily manageable. 
In three spatial dimensions, solutions to equation l|tjtj|) include ellipsoids, disks, biconcave and torus-like (or "wheels" 
in the terminology of @) forms 0,@- We now wish to consider if one can use these toroidal fluid configuration to 
model strong gravity configurations with the same topology. In particular, we consider a torus with a large radius R 
and small thickness r + (the reason being that it is in this regime that one really expects a "pure" torus). 

We will first look for torus- like solutions which start smoothly (i.e. <f» = 0) at some distance R\. We will call 
these "type A" toroidal solutions. The outer radius will be denoted R2, where <fr = — 00. In this case, we can easily 
determine the constants po , C and write Ij66|l as 



-D-l ud-i 



r 1 - 1 P n 2 (w D+1 - R? +1 ^ 



Pn 



,D-2, 



T(D - 1) 



where 



oD-l pD-1 

Po A = R D - 2 



2T(D + 1) y/TT^' 

P n 2 - 



T(D-l) " 2 2T(D + 1) 

We have numerically integrated equation Q79JI for several values of the parameters, and we find the law 

R 2 (R2- Rif {D + 1) A ~ 3.57 , 



which is independent of dimensionality D. Here 



A 



P n 2 

2T(D + 1) 



(78) 



(79) 



(80) 



(81) 
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FIG. 5: Typical equilibrium figures of rotating drops. The (blue) dashed line corresponds to a spherical surface with Ef, = 0. 
The (magenta) dotted line corresponds to Ef, = 2 and the (black) solid line corresponds to the limit E;, = 2.329. In all three 
cases the equatorial radius was chosen to be unity. 




0.00 0.02 0.04 0.06 0.08 



FIG. 6: Typical "smooth" Type A torus-like equilibrium shape. We have taken D = 4 and A = 3.5, Ri = 99.9549, R 2 = 100. 



The meridional cross-section for a typical type A torus is displayed in Fig. As shown by Gulliver [{J and Smith 
and Ross |45j there are several possible families of toroidal fluids (some of which have been experimentally observed 
already, see 3]). We have just described one of those families, that starts smoothly at R±. 

We can find a second family, the type B toroidal shape, which is more similar to the gravitational Dyson rings 
which will be described in Section^ What we now look for is a family which starts "abruptly" at R\. Our numerical 
experiments seem to suggest a very rich behavior for the members of this family. For very small A, we find the 
following approximate scaling law, 

1 QR 

r + R 2 A~—— , R=R 2 ~R X . (82) 
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and the meridional cross-section of a typical member of this family is shown in Figure Such a family seems to 
have been discussed already by Appell [46|. Let's pause a moment on this scaling law. If, following the suggestion in 




I / / . 1 

0.0 0.2 0.4 0.6 0.8 20000.00 20000.01 



FIG. 7: Typical Type B torus-like equilibrium shape. We have taken D = 4 and A = 1CT 7 , Ri = 2x 10 4 , R 2 = 2.00000099 x 10 4 
We chose — 3*- 



T(D-l) 



|22|. we try to model gravitating black bodies with fluids held together by surface tension, the first law suggests [22] 
that the surface tension is the equivalent of the Hawking temperature (up to a multiplicative constant). Now, the 
temperature of a black hole scales as 1/r, where r is some typical radius. For large "black rings" the relevant radius 
is the thickness r + of the ring. Thus we set T ~ l/r+. Using l|81[l and (|82|l we get 

M 

^ 2 ^3-l 5 (83) 

where we replaced p = M/V. 

On the other hand, the mass per unit length of such a "black torus" should be given by the mass of a D — 1 
Tangherlini black hole J37J (i.e., a higher dimensional non-rotating black hole): 

{D - 2)C D GM ~ 4tt(L> - 2)C D - 2 Rr+~ 3 . (84) 

Using the "constitutive" relation l|84|) for black objects and J ~ Q.M R we get finally 

J 2 



APR 2 



1 . (85) 



Fortunately, there is an exact solution of Einstein's equations describing a general relativistic black ring, in five 
dimensions (D = 4) [4?il |: the black ring solution by Emparan and Reall JiJ . The black ring solution is described by 
two parameters v and R. They are characterized by having mass and angular momentum equal to [4oL l4l|: 

4G L l + v y ' 

ttR 3 V\i>(1 + A) 5 / 2 
J 2G~ L (l + v) 2 ■ (87) 

Here, the parameter A = 2^/(1 + v 2 ) for black rings. Large rings can be obtained by setting v — > 0, in which case R 
can be looked at as the radius of the ring. For very small v we then have 

J V2 



MR 3 ' 

This is consistent with relation (|85|l obtained in the model of a fluid held together by surface tension. The case D = 3 
seems to be special for in three dimensions there are no black strings (all black holes must be topological spherical in 
this case), whereas the surface tension model seems to allow for rings; 
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VII. STABILITY ANALYSIS OF A ROTATING DROP IN CASE I. ONE NON- VANISHING ANGULAR 

MOMENTUM 



To complete our discussion on drops held by surface tension, we consider the generalization to D-spatial dimensions 
of Figure ^ the bifurcation and stability diagram. Here we restrict ourselves to case I, i.e. only one non- vanishing 
angular momentum, corresponding to rotation in the plane x 1 — x 2 , and we will necessarily leave many topics for 
future research. For instance, all asymmetric figures will be left out of this exercise. We will follow the approach of 
Chandrasekhar |5j], based on the formalism of virial tensors. The special case in which the rotation is zero, i.e., a 
spherical drop, is considered in Appendix^where we compute the oscillation modes of a Z?-dimensional incompressible 
fluid sphere held by surface tension. 

We have found the basic equations describing the equilibrium shape of rotating drops in the previous Sections. We 
would now like to investigate if these shapes are stable. The formalism to do this was laid down in Section II V CI and 
equations (|55|I - H58|) are the most general set of equations describing perturbed rotating drops. 

From now on we will restrict our analysis to a special kind of perturbations, termed "toroidal" modes by Chan- 
drasekhar [5J. The reason is that these modes are expected to play a role in the bifurcation and stability analysis. 
Future work should include an analysis of the other modes (termed "transverse-shear" and "pulsation modes"). We 
refer the reader to the work by Chandrasekhar p| for further details on the formalism. The toroidal modes are defined 
by 



= ax 1 + (3x 2 
£ 2 = fix 1 - ax 2 
Z A = 



(89) 



with a, (3 constant parameters. 

As shown in the Appendix El if a perturbation has the form (|89|) . the only non- vanishing virial integrals are 



r i;i 



-^2:2 



= irpa 
npf3 



G'd-2 I 3r,, m-D-2 , 

w [j[vj)\ aw 



D-2J Q 

Cd-2 



D 



I" w 3 [f{vj)] D ~ 2 dm 
Jo 



(90) 
(91) 



where vo — a is the intersection of the surface of the drop with the x 1 — x 2 plane. These quantities satisfy the equations 
(|55|l . which have the form 



and, combined, become 



lu 2 t 12 



w 2 r 1;1 - 2w«r 2;1 


= <56 n - 


Hf2 2 r n + 


an, 


lu 2 T 2 '< 2 + 2luQT U2 


= se 22 - 


i- n 2 r 22 + 


6n, 


w 2 r 1;2 - 2cj^r 2;1 


= se 12 - 


f n 2 r 12 , 




w 2 r 2;1 + 2luQt U2 


= se 21 - 


f n 2 r 21 




- r 22 ) - 2luQt 12 


= se 11 - 


se 22 + n 


2 (r n 




= 2se 12 - 


f 2ft 2 r 12 . 





(92) 
(93) 
(94) 
(95) 



(96) 
(97) 



In order to solve equations H96|) . 1971) . where the virial integrals are given by (|9U|I . (|91|l . we have to compute the surface 
integrals <56 12 and S6 11 - 5& 22 . 

In the case of perturbations of the form they are (see Appendix IDl) 



12 



se 



se 11 - se 22 



T 



(3nC 



D-2 



dw^w^Qivj) [/(ro)] 



D-3 



-TanCD-2 / dw4>w Q(w)[f(vj)\ 
Jo 



D-3 



(98) 
(99) 



where 



4>'{4> 2 -2) 

W 2 V(l +0 2 ) 5 / 2 



+ . 



2 + 30 2 



D-3 



7(1 + 2 )3/ 2 



Z (1 + 2 ) 3 / 5 



(100) 
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Therefore equations (|96|l . (|97() are equivalent to 

Lu 2 a - 

uj 2 (3 + 2Ma = -23P + 2fl 2 p (101) 



u 2 a-2ujnp = -23a + 20 2 a 



where 



„ _ T( D - 2) J a 0CT 4 Q(tu)[/(tu)] g - 3 rf n7 
Setting = — a 2 , the algebraic system (|101fl admits solutions 



cr = ft ± ^23 - ft 2 . (103) 

A neutral mode of oscillation (<r = 0) occurs at f2 2 = J, but (in the absence of dissipative mechanisms as stressed in 
@) the drop is still stable in the vicinity of this point. If 

Q 2 <23, (104) 

then a is real and the perturbation is stable; if the rotation frequency of the drop is larger, then the drop is unstable 
against toroidal perturbations. To ascertain stability or instability, we thus have to use the results of section TV Al for 
the equilibrium figures of rotating drops and insert the numerically generated shapes z = z{vo) in (|102|l and H103[l . 

If we measure 137 and z in units of a and fl in units of w we get 



0= v/s;± v/2J-E a , (105) 

where the dimensionless quantity J is 

= QD - 2) £ ^Q{w)[f{w)] D ^dw 

32 j>3 [/rf - 2(fo ■ ^ 

The results of solving (|105fl are displayed in Fig. [8] for D = 4, which is a first step towards reconstructing the whole 
bifurcation diagram. In Fig. [SI o\ is (for real solutions) the largest of the two roots in (|1U5|I , i.e., as long as 2 J > E a . 
For 2 J < S a , we display only one root (the other root is the complex conjugate of this) by plotting the real and 
imaginary components of a. 

In three dimensions (D = 3) it was found by Chandrasekhar that rotating drops exhibit a neutral mode of 
oscillation at £ = 0.4587. He also found that rotating drops become unstable for S > 0.84. We find that for general 
D the results are similar. For instance, in D = 4 there is a neutral mode of oscillation at E a = 0.46, and above the 
critical E a = 0.86 the drops are unstable. The shapes of the drops at these values of S are displayed in Figure 0] 
Results for higher D are very similar. Our numerical method only allowed to probe dimensions up to D = 6, for 
larger D we loose accuracy and the results cannot be trusted. 



VIII. STABILITY OF FLUID CYLINDERS HELD BY SURFACE TENSION: THE 

RAYLEIGH-PLATEAU INSTABILITY 

To end our discussion on fluids held by surface tension and also as a warm up exercise for the up-coming Sections, 
we will review the results of |2^| , describing a simple capillary instability - the Rayleigh-Plateau instability |47l Esj - 
as a good model for the Gregory-Laflamme instability. Most calculations were omitted in [23| so we present the full 
detailed computations here. 



A. Rayleigh-Plateau instability in higher dimensions 

We take an infinite non-rotating cylinder of radius Ro held together by surface tension (therefore without self- 
gravity). Neglecting the external pressure, the unperturbed configuration is worked out easily to yield an internal 
pressure Pq = (D~2)T/Rq, where D is the total number of spatial dimensions. We perturb slightly this configuration, 
according to 

r = R. + ee ult+lkz . (107) 
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FIG. 8: The dependence on the quantity E of the characteristic frequencies of the toroidal modes of a rotating four-dimensional 
drop in D = 4. The diagram is very similar to the three dimensional one (D = 3) presented in Figure 3 of Chandrasekhar 
As a check on our numerical code, we verify that in the non-rotating limit our numerical results agree well with the analytical 
prediction a = \[D + 1/2. 



The equations governing small perturbations are 

Po + SP = Tdivn, (108) 

= -gradT, (109) 

divu = 0, (110) 
SP 

Y = — 7 (HI) 
P 

where Pq is the pressure in the undisturbed cylinder and u is the fluid velocity. For instance, u r is the radial component 
of the velocity vector. Equation (|109|l yields 

u r = --j-T, (112) 

uj dr 

u z = — T. (113) 

Equation 1110|l implies that 

du r D — 2 du z 

- f L + u r + —^ = 0. (114) 

ar r dz 

Substituting in this the relations l|112Jl and 1)113(1 we get 

d 2 T D-2dT ,„ , , 

— + — -fc 2 T = 0, 115 

dr z r dr 

with (well-behaved at the origin) solution 

T = Cr 3 ^- D / 2 I a {kr), (116) 

where a = (D — 3)/2. By relation l|112|) we get 

u r = -C l——I a (kr) + kl' a {kr) \ , (117) 

where I' a (x) stands for the derivative of the Bessel function with respect to the argument x. Now, for this velocity to 
be compatible with the radial perturbation (|107J) one must clearly have 

ee ^t +tkz = -C^L_ (±__ Ia{kRo) + kI' a (kR )\ . (118) 
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With the displacement (|107f) we have to first order in e 



— = — £- (D - 2 - (fci? ) 2 ) ee" t+lkz . 
P RqP 



Using (fTTHf) and (|TT%|) together with fTTTjl . (fTTSf) we finally find 



2, ,2 



r(w 7 «( fci? o) + fc^(fc^o) 



/0 



(U - 2 - (fei? ) 



(119) 



(120) 



where I' a (x) stands for the derivative of the Bessel function with respect to the argument x. This is the sought relation 
yielding the characteristic frequencies of an infinite cylinder of fluid held by surface tension. Formula l|120|) constitutes 
a generalization of a classic result by Rayleigh ^t|> an d it predicts that for kRo < y/D — 2 the cylinder is unstable 
against this kind of "beaded" perturbations. In Figure we show ui as function of wavenumber k for several space 
dimensions D (actually what is shown in the plots are dimensionless quantities loRq and fci?o)- In Section llXI we will 




FIG. 9: The Rayleigh-Plateau instability of a hyper-cylinder held together by surface tension, in several dimensions. Here the 
effective density were chosen to match those of a higher-dimensional non-rotating black hole. The instability gets stronger as 
the spacetime dimensionality increases, and the critical wavenumber (for which uo = 0) increases with D. 

discuss this instability in more detail, when we compare it to its self- gravitating counter-part. 



B. Rayleigh-Plateau and Gregory-Laflamme 



In General Relativity, there is an interesting class of black objects (objects with an event horizon) without spherical 
topology, the black branes 39] . A broad class of these objects are unstable against gravity, in a mechanism known 
as Gregory-Laflamme (GL) instability |44j . This is a beaded instability that makes any small perturbation with 
wavelength A of the order of, or larger than, the radius of the cylinder Rq grow exponentially with time. It is a very 
robust long-wavelength instability. In |22j it was proposed to mimic the GL instability with the Rayleigh-Plateau 
instability, based on the suggestive appearance of the second law of black hole mechanics which endows the horizon 
with an effective surface tension. It was found that all aspects of the GL instability could be accounted for by the 
Rayleigh-Plateau instability in fluid cylinders. In some sense, the effective surface tension can mimic all gravitational 
aspects of black strings. If this is indeed correct, we should be able to predict the effects of rotation and charge on 
the GL instability, as discussed in the following. 



1. Effects of magnetic fields 

Even though not pursued here, a generalization for charged and/or rotating cylinders is of interest. In any case, we 
can lean on the four-dimensional results (D — 3 spatial dimensions), for which the charged and rotating cases have 
been worked out [43,|^0|. A magnetic field has a stabilizing effect j4J|, which is in agreement with the fully relativistic 
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results for the Gregory-Laflamme instability, where it has been shown that adding charge to the black strings stabilizes 
it 01 ■ There is good evidence that extremal charged black strings are marginally stable, and indeed in the capillary 
case (cylinder held by surface tension) it is possible to completely stabilize the cylinder if the magnetic field has 
a certain strength (49j (dependent on the conductivity of the fluid), so this may simulate the marginal stability of 
extremal black strings very well |44|. Nevertheless an extension of the calculations presented here to the case of 
magnetically charged cylinders with surface tension is needed. 

2. Effects of rotation 

The effects of rotation on the capillary Raylcigh-Plateau instability can be discussed more intuitively. The cen- 
trifugal force, scaling with r, contributes with a de-stabilizing effect since there is an increased pressure under a crest 
but a reduced pressure under a trough. So in principle, the threshold wavelength should be smaller |Hcj . i.e., the 
instability should get stronger when rotation is added. It is natural to expect that rotation will lift the degeneracy 
of the modes and that unstable modes will in general be non axi-symmetric. What can we expect for the general 
relativistic GL instability once rotation is added? A very naive reasoning is that since extremal charged black strings 
are marginally stable then perhaps the same happens with rotation, and so rotation has a stabilizing effect. This is 
incorrect as a simple argument shows. Consider a rotating black string in D + 1 spacctime dimensions (with again D 
being the number of spatial dimensions) of the form 

.9 = Kcrr D x R , (121) 

where Kerr 15 is the metric of higher dimensional Myers-Perry rotating black hole in D spacetime dimensions, with 
only one angular momentum parameter (case I) . If we focus on D > 4 the angular momentum of these black holes is 
unbounded and one can show that for very large angular momentum 

Kerr 15 ~ S^~ 2 x R 2 , for large rotation (122) 

where S 1 ^" 2 is the metric of a higher dimensional Schwarzschild black hole, often referred to as Tangherlini |3?1 | black 
holes. We thus finally get 

g = S°- 2 x R 3 , (123) 

This is the metric of a non-rotating black membrane in D + 1 spacetime dimensions extended along the 3 dimensions. 
But this is precisely the geometry considered by Gregory and Laflamme 44], and so it is unstable. While this 
argument says little about the strength of the instability as function of rotation, it does predict that rotating branes 
are still unstable. Putting this and the fluid membrane argument together we predict that rotation makes the GL 
instability stronger. One can also use an entropy argument, similar to the one in to show that rotating black 
strings should be unstable. Again, an extension of the calculations presented in Section lYlII Al to the case of rotating 
cylinders with surface tension is needed. The D = 3 case, studied in |53| . indicates that other effects may come into 
play with rotation. For instance, other instabilities seem to manifest themselves for high enough rotation, so this is a 
particularly exciting topic to consider. 

Part II 

Fluids held by (Newtonian) self-gravity 

IX. STABILITY OF SELF-GRAVITATING FLUID CYLINDERS: THE 
DYSON-CHANDRASEKHAR-FERMI INSTABILITY 

It is only appropriate that we start the study of self-gravitating fluids by examining the (self-gravitating) counterpart 
of the cylinder studied in the last Section lYlIII 

A simple consistent solution can be found in £>-dimensions, which corresponds to a static self-gravitating infinite 
cylinder made of fluid with a density p. For the moment we will not consider rotational effects (these will be discussed 
later), and we will focus on incompressible fluids. The stability of such a solution for the case of 3 spatial dimensions 
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was studied for the first time by Chandrasekhar and Fermi [5l| (see also |49j), who showed that such a solution 
is unstable against cylindrically symmetric perturbations (those for which the cross-section is still a circle, with a 
z— dependent radius), and stable against all others. If the radius of the cylinder is Ro, this "varicose" instability (in 
Chandrasekhar's [4!| terminology) sets in for wavelengths A > 27ri?/1.0668. 

A similar phenomenon had however been studied already in 1892 by Dyson |18|. who investigated the stability 
of self-gravitating rings of fluid. In his terminology, the ring is stable against fluted (the ring remains symmetrical 
about its axis but the cross-section is deformed) and twisted (the cross-section remains circular but the circular axis 
of the ring is deformed) perturbations. He also showed that sufficiently thin rings are unstable against a third kind of 
perturbation, so-called beaded perturbations (those for which the circular axis is undisturbed, but the cross-section 
is a circle of variable radius). But these correspond precisely to Chandrasekhar-Fermi type of perturbation (with the 
added bonus that Dyson effectively shows that the instability goes away if the length of the cylinder or ring is less than 
the critical wavelength), and thus we shall refer to this classical fluid instability as the Dyson-Chandrasekhar- Fermi 
(DCF) instability. 



A. A self-gravitating, non-rotating D-dimensional fluid cylinder 

We will now generalize the DCF instability for an arbitrary number of dimensions. Consider an infinite cylinder 
of incompressible fluid, with density p and radius R. To fix conventions we consider spacetime with 1 (time)+l (axis 
of cylindcr)+Z? — 1 = D + 1 total number of spacetime dimensions. In the following the coordinates r , z have their 
usual meaning in polar coordinates. The gravitational potential *6(x) inside and exterior to the cylinder satisfies 
respectively 

V a fl(x) 4 = -(D-2)C D Gp, (124) 
V 2 B(x) e = 0. (125) 

Now, in D spatial dimensions with cylindrical symmetry we have 



v 2 = J__±( r D-2±\ + ^_ (126) 

r D 2 dr \ dr J dz 2 



Since we consider a cylinder with constant radius (along the axis in the z-direction), the last term drops off and we 
get the following solution (we demand that both the potential and its derivative be continuous throughout) 

B(x)$ = - ^V , (127) 

B(x) e = - {D ~^° D GpRl logr/Ro - E^lc D GpI% , D = 3. (128) 

_ (D 2)C D G P R 2 ( _D-1_ 2<- 3 \ 
° Wo _ 2(D-1) l(D-3) (D-3)r D -* r [ ' 



Hydrostatic equilibrium demands that the pressure Pq satisfies 

dr \ p 

which yields, upon requiring Pq{t — Rq) = 0, 



4- f^r -B(x) | = » CI.".") 



B. Axisymmetric perturbations 

Suppose now the system is subjected to a perturbation (see Fig. I10[) such that the boundary is described by 

r = i? + ee ut+lkz . (132) 
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\=2%/k 

FIG. 10: Black strings and fluid cylinders are unstable to perturbations on the extended dimension, i.e., along the axis of 
the cylinder. Ripples propagating along this axis grow exponentially with time for wavelengths of order of the radius of the 
cylinder. 



The equations governing small perturbations are 

(133) 

radT , (134) 

(135) 

- SB{xY . (136) 

Here u stands for the velocity of the perturbed fluid. Let's solve for the first of these equations. When solving for the 
gravitational potential, the last term in (|126[) cannot be neglected. We find 

£(x) e = B{x)l + Aee" t+lkz r z t 2 - D ' 2 K a {kr), (137) 
£(x) 1 = B(x)j + Bee^ t+lkz r^ 2 - D / 2 I a (kr), (138) 



V 2 <5£(x) 


= o, 


du 




at 




divu 


= o, 




SP 


T 






P 



where 



D — 3 

a = —-±, (139) 



and K a and I a are Bessel functions of argument a. The constants A and B are determined by requiring continuity 
of the potential and its derivative on the boundary. By keeping only dominant terms in e we find 

A = {D-2)C D GpR I a (kR Q ), (140) 
B = (D -2)C D GpR K a (kR ) . (141) 

Therefore we have 

5B{*) 1 = e(D - 2)C D GpR K a (kR )e UJt+ikz r 3 / 2 - D / 2 I a (kr) . (142) 
We also have u ~ e wt + tkz ; therefore Eq. I|134|) yields 

u r - — — y-T , (143) 
uj ar 

u z = — T. (144) 

(Here u r , for instance, is the radial component of the velocity vector). Equation (|135[l implies that 

du r D-2 du z 

— 1 u r H — j— = . (145) 

ar r dz 

Substituting the relations (|143f) and (|144f> we get 

d 2 T D-2 dT 



dr 2 r dr 



fc^T = 0, (146) 
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with (well-behaved at the origin) solution 



By relation i|143|) we get 



T = Cr^ 2 - D / 2 I a {kr). 



(147) 



3/2-B/2 / 3 _ D 

n r = -C [ — I a {kr) + kl' a (kr) 



2r 



(148) 



where I' a (x) stands for the derivative of the Bessel function with respect to the argument x. Now, for this velocity to 
be compatible with the radial perturbation (|132f> one must clearly have 



p 3/2-D/2 , „ _ n 
ut+ikz _ _nf}^_ 1 '' 1 

' 2Rn 



I a (kR ) + kI' a (kR ) 



(149) 



Our problem is completely specified once we impose the boundary condition that the pressure must vanish at the 
boundary, 



an ap 

SP + e^- =0,^ pT + pSBixy + e^p- = 0. 
or or 



This gives us finally 



^ = (3-^(^o) + ^(^o) ( Ia{kRo)Ka{kRo) _ _J_ 

(D-2)C D G P 2I a (kR ) V 0J ay aj D-l 



(150) 



(151) 



In Fig. ^]we plot the dispersion relation oj(k) for several values of the dimensionality D, using Eq. I|151fl . When 




FIG. 11: The DCF instability of a self-gravitating hyper-cylinder in several dimensions. Here the effective densities were chosen 
to match those of a higher-dimensional non-rotating black hole. The instability gets stronger as the spacetime dimensionality 
increases, and the critical wavenumber (for which uj = 0) increases with D. 

D increases the maximum ui grows which means that the instability gets stronger with space dimensionality. The 
threshold wavenumber k c , defined as the wavenumber for which w = 0, also increases with D. In Table [i] we present 
the threshold wavenumber for selected values of D. 

For large D we can obtain an analytical expression for the threshold wavenumber. In fact we can expand the Bessel 
functions as 1551: 



I a (kRo) ~ (fci? /2)« ( 1 + {kRo)2 



a\ 4(a + l)! 
K a (kR ) ~ l W 2)-°((a-l)!- ^- 2 ^ 2 



(152) 
(153) 
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TABLE I: Dimensionless threshold wavenumber kRo for the DCF instability of a higher dimensional fluid cylinder, and the 
corresponding threshold wavenumber for Rayleigh-Plateau and Gregory-Laflamme instability (these last data are taken from 
[H^H . The cylinder has D — 1 dimensions with radius Ro plus an extended direction along z. For D = 3, not shown in the 
Table, we get kRo = 1.06672, in complete agreement with Chandrasekhar and Fermi's |51j result. 



D 


4 


5 


6 




7 


8 


9 


49 


99 


DCF 


1.41 


1.70 


1.96 




2.20 


2.41 


2.61 


6.85 


9.85 


Rayleigh-Plateau 


1.41 


1.73 


2.00 




2.24 


2.45 


2.66 


6.78 


9.80 


Gregory-Laflamme 


0.876 


1.27 


1.58 




1.85 


2.09 


2.30 


6.72 


9.75 


We get 


K a {kR Q )I a (kR ) - 


1 

D - 1 


2 

' 


(kR ) 2 
4q 3 


4> kR — 


Vd. 




(154) 



A comparison between these results and the ones in Section lVlIIl is instructive. The Figs. I§1 and ITTI show a remarkably 
similar behavior for both these instabilities. The threshold wavenumber (for which u = 0), shown in Table [i] is also 
very similar to the corresponding one in both the Gregory-Laflamme and the Rayleigh-Plateau instability. Finally the 
D — > oo limit of the three instabilities is the same, kRo ~ V~D, in the sense that the threshold wavenumber coincides 
in this limit. 



C. DCF and Gregory-Laflamme 

In [23, the authors proposed to mimic the GL instability with the Rayleigh-Plateau instability, based on the 
suggestive appearance of the second law of black hole mechanics which endows the horizon with an effective surface 
tension. What we now propose is that, while the fluid membrane analogy is extremely powerful (and perhaps more 
trustworthy when it comes to mimic black objects), the Newtonian classical counterpart of the GL instability is 
the DCF instability. The GL instability should be the limit of the DCF instability when the density is so large 
that the cylinder is in fact a black string. Not only are the details of the DCF instability very similar to the GL 
instability (for instance, compare Figure ITTI with Figure 1 of 0]) but for large D we get exact agreement for the 
threshold wavenumber. In Tabic [i] we compare the threshold wavenumber for the three different types of instabilities 
(DCF, Rayleigh-Plateau and GL). The agreement is indeed very good, thus providing a convincing evidence for our 
conjectured analogy between DCF/Rayleigh-Plateau instabilities for fluids and GL instability for black strings. 

It has been shown by Chandrasekhar and Fermi [EJ (see also Chandrasekhar ^^|) that the DCF instability is 
rooted in the sign change of potential energy for these varicose or beaded perturbations. But this is precisely what 
happens for the capillary of Rayleigh-Plateau instability (see Section IVlIII and discussion in [23 ), which explains why 
both of them yield so similar results. This is another pleasant indication that these instabilities are related. Based on 
this agreement we can discuss what to expect for rotating and charged cylinders and discuss the significance of the 
results for the GL instability. 

D. Effects of rotation and magnetic fields 

Even though not pursued here, a generalization for charged and/or rotating cylinders is of great interest. In any 
case, we can lean on the four-dimensional results (D — 3 spatial dimensions), for which the charged p^, l5ll l5t | and 
rotating cases [57|, |5l| have been worked out. 

It has been shown that a magnetic field has a stabilizing effect, reducing the threshold wavenumber, i.e., increasing 
the minimum wavelength at which instability occurs. This is in agreement with the fully relativistic results for the 
Gregory-Laflamme instability, where it has been shown that adding charge to the black string stabilizes it . There 
is good evidence that the extremal case is marginally stable, and this could be perhaps the only real signature of 
an event horizon |78|. because in the horizonless, Newtonian case, "no magnetic field, however strong, can stabilize 
the cylinder for disturbances of all wavelengths: the gravitational instability of long waves will persist" |4^. As we 
discussed in Section IVIIII a capillary model may be more accurate to describe extremal objects. 

Focusing now on the effect of rotation, it has been shown [EtJ (in three spatial dimensions) that rotation has a 
destabilizing effect. This agrees with the phenomenological arguments in [23] for cylinders held by surface tension 
(see also |5(|) of the capillary model. It also agrees with the arguments in Section IVIIII suggesting that rotating black 
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strings are unstable and thus a solid conjecture is that the Gregory-Laflamme instability of black strings gets stronger 
if one adds rotation. Stronger here means the threshold wavenumber decreases (the threshold wavelength increases). 
An interesting extension of the results presented here would be to consider a Z?-dimensional rotating cylinder of fluid. 



X. SELF-GRAVITATING TORUS AND BLACK RINGS 



Toroidal shapes are thought to be common configurations in the universe. In some of them these rings or tori are 
associated with a massive central object, such as Saturn's rings, whereas in others the central object plays almost no 
role at all (see and references therein). Toroidal configurations in full general relativity are also important (see 
|62j|V and they sometimes appear in the form of temporary black toroids, having an event horizon, in simulations of 
collapse of rotating matter [&| . 

The problem of a self-gravitating torus in Newtonian theory was investigated by many authors, the first studies 
belonging to Poincare, Kowalewsky E2 and Dyson Dyson E3 considers a self-gravitating ring of radius R and 
thickness 2r + (i.e. the cross-section of the ring is a circle with radius r + ). For small r + /R he finds a power series 
solution (the leading term had already been derived by Poincare and Kowalewsky |Tlj ) to the equilibrium problem as 



nG L p 



£N"H)-£M"B)*° 



(r+/Rf : D = 3. (155) 



To generalize this result to arbitrary dimensions we first observe that the gravitational potential of a very large ring 
of radius R can be obtained by using as gaussian surface a torus. We get, for D > 3, 



(D - 2)C D GM 
Att{D - 3)C D - 2 Rri 



g(*)~- , ,r :v 3 > ^ 



where we expressed everything in terms of the physically more significant quantity M, the total mass of the ring. For 
rings with a large radius R the hydrostatic equilibrium condition can be written as HI Eg tl 2 R 2 ~ B(x). This gives 



us 



n 2 R 2 (D-2)C D M 

4n(D - 3)C D -2Rr°-* 



It seems plausible that this formula also describes general relativistic objects, as long as one uses an "equation of 
state" to deduce r + in terms of the other quantities. This is because we are considering large rings, without really 
appealing to the internal structure of the object. For large rings, the dominating force determining equilibrium 
should be well-described by Newtonian forces, even though general relativistic effects become important "locally" and 
therefore it is the latter that determine r + (R, M). To plug-in the information that we are dealing with a black object, 
we use the relation between r+, M and R. For a black ring, the mass per length should be given, to a good accuracy, 
by the mass of a D — 1 Tanglier lini black hole [37| : 

{D - 2)C D GM ~ 4tt(L> - 2)C D - 2 Rr+~ 3 . (158) 

Using l|157|) and re-expressing everything in terms of the angular momentum J ~ M R 2 fl we get 

J 2 



M 2 R 2 



(159) 



We can compare this prediction against the exact general relativistic black ring solution |40| in five spacetime dimen- 
sions (D = 4). For large black rings, Eq. JSSJ), jjr = 3 > m verv good agreement with the Newtonian prediction 
(|159fl . Moreover these arguments also lead to conjecture the existence of black rings in higher dimensions. Again 
the case D = 3 seems special for there are Dyson rings in this case (actually that's the case where they were found), 
albeit with a peculiar logarithmic term. As we remarked at the end of Section IVII there are no black rings in D = 3 
so this is an instance where the analogy breaks down. To end this section, we note that another accurate and simple 
model was put forward by Hovdebo and Myers |64| to describe black rings in quasi-Newtonian terms, starting from 
the black string solution. Their conclusions are basically the same as ours, with the exception that instead of working 
with Newtonian forces, they extremize the entropy of the solution. 
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XI. SELF-GRAVITATING, ROTATING SPHEROIDS: THE MACLAURIN SEQUENCE 

To complete our discussion of self-gravitating objects, we consider the generalization to D-spatial dimensions of the 
simplest rotating objects: the MacLaurin sequence. Here we restrict ourselves to case of one non-vanishing angular 
momentum (case I in the notation of Section^), corresponding to rotation in the plane x 1 — x 2 ; the coordinates are 
then decomposed as (J7J 

x l = {x a ,x A ); a =1,2; A = 3,...,D. (160) 

We will follow the approach of Chandrasekhar ,12), based on the formalism of virial tensors. The special case in which 
the rotation is zero (i.e., a self-gravitating globe), is considered in Appendix 151 where we compute the oscillation 
modes of a D-dimensional incompressible fluid sphere. 



A. The MacLaurin sequence in three dimensions 

Here we briefly review the classical calculation of the MacLaurin sequence due to Chandrasekhar. We refer to [l^ 
for further details. 

Let us consider a matter distribution, inside the volume V, with density p(x). Its gravitational potential is 

B{x) = G ( T^Xr^x' . (161) 
Jv \x-x'\ 



We define the tensor 



which satisfies Bu — B. We also consider the potential energy 



B, , ( .7 > = (V / p(5>) {X% f" ){X ' X ' J) d a x' , (162) 



and the tensor 



The following relation holds: 



Indeed, 



JB=~ I f >\.7)Bir ur.r. (!(,:-!) 



Mij = ~\j Pi s ) B H W d * x ■ ( 164 ) 



IBij = J^p^-x^d 3 x. (165) 



p^x^d 3 x = -G [ d 3 x [ d?x'p(x)p(x') X ^ i* xi = G [ d 3 x' [ d 3 xp(x')p(x) X '\ x'i 
v dx l Jy Jy Fy '\x-x'\ 3 J v J v ; |f-f'| 3 

= ~G I d 3 x f d 3 x'p{x)p^) {xl -^ ){x '- X ' ]) . (166) 
* Jv Jv I a; — x i 

A consequence of this formula is that lBij is symmetric in its indices. 
We then have the equations of motion 

BP BE 
dP dB 

z = -^t + Pjt, 167 
oz oz 

where z = x A — x 3 . 
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At equilibrium, 



n2 a dP dB 



-T + P 9 /=0- (168) 
az az 

We compute the virial integrals by multiplying equations (|168[) by x J and integrating in the volume domain V . The 
resulting integrals are: 



O 2 f d 3 xpx a x j = n 2 i aj , 
Jv 

f> x ^ x i = 5 ij f d 3 xP = SiiIl. 



v 



dx l 



Thus we get virial equations 



8B 

d 3 xp—x : > = JBij. (169) 

fl 2 Iab + IBab = —HSab , 

n 2 i a3 + m a3 = 0, 

m 3a = (therefore IB a3 = hs = 0) , 

IB 33 = -II. (170) 
The inertia tensor at equilibrium is then, in general, 

/ hi h2 

hi hi o ! uni 
\ o 0/ 3 3 



hj 



We note that equilibrium does not require necessarily hi = Iii- 
From the virial equations we get 

JBxx + r» 2 /n = IB22 + ^ 2 / 22 = IB 33 ■ (172) 

There are two possible solutions of equations l|172|) . The first is the simplest one, with 

In = I22 , 2B U = IB 22 , (173) 

namely, with axial symmetry; the rotation frequency is given by 

n2= IB 33 -IBn 
hi 

and can be easily computed; this solution corresponds to the MacLaurin sequence of oblate axially symmetric ellipsoids. 
The second, more complicate, has hi 7^ I22, and corresponds to the Jacobi sequence of tri-axial ellipsoids; this solution 
exists only if the total angular momentum exceeds a critical value. 



B. The MacLaurin sequence in higher dimensions 

The generalization of the three dimensional formalism to higher dimensions is more or less straighforward. Let 
us consider a matter distribution in D dimensions, inside the volume V, with density p(x). Again, we define the 
gravitational potential 

B(x) = G [ ^ Pi fJ n _ 2 d D x' . (175) 

and 



J iJ 



(*) es f P(af) ( - 0§^ X ' J) d»x> , (176) 
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which satisfies Bu = B. We also consider the potential energy 

m=-\ { P {x)B{x)d D x. 
2 Jv 

and the tensor 

IB ij = -~Jp(x)B ij (x)d D x, 
We can check that the following relation, which is the extension of (|165|l . holds: 



p-—^x J a x. 

v ox 1 



(177) 



(178) 



(179) 



Indeed, 



dB 



v dx i 
D-2 



p^x j d D x = -{D-2)G / d D x / d D x'p{x)p(x')-^ — ^x 1 = (D - 2)G / d D x' / d D xp{x')p{x) 



x — X 



\x — x 



'IID- 



G I d D x d D x'p(x)p(x') 



, Ax 1 -x' i )(x j -x' j ) 



\x — X 



i\D 



(180) 



Thus IBij is symmetric in its indices, as it was in three dimensions. 
The equations of motion yield 



^ + n 2 px a + 2 P ne ab x b +p^-, 
ox a ox a 

dP dB 



At equilibrium we get, 



dx A dx A 



il 2 x a - — + — 

^ dx a ^ dx a 

dP dB 



(181) 



dx A dx A 



(182) 



We compute the virial integrals by multiplying equations (|182f> by x 3 and integrating in the volume domain V. The 
resulting integrals are: 

fl 2 I d D xpx a x J = Q 2 I n 



v 



dP 

dx l 



- I d D x^-xJ = S ij I d D xP = 5,~Il. 



I v dDxp S xj = ( £> - 2 ) jb * 



(183) 



Thus we get the virial equations 

n 2 i ab + [d - 2)iB ab = -m ab , 

Q 2 I aA + {D - 2)IB aA = 0, 

(D ~ 2)IB Aa = (therefore IB aA = I aA = 0) , 

(d-2)ib ab = -m AB . 

The inertia tensor at equilibrium is then, in general, 

/ hi h2 \ 

J33 



(184) 



V 



(185) 
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We note that equilibrium does not require necessarily In — I 2 2- 
From the virial equations we get 

(D - 2)2B U + n 2 I n = (D- 2)5? 22 + fl 2 I 22 = (D - 2)1B 33 



(D - 2)IB DD . 



(186) 



The simplest solution of Eqs. (|186fl is the one with axial symmetry: 



hi — I22 , 



IB n = IB22 



(187) 



which is the generalization of MacLaurin sequence; the rotation frequency is given by 



VI 2 = (D - 2) 



hi 



(188) 



and can be easily computed. But equations l|18tj|) leave open the possibility to have other, non axisymmetric equi- 
librium configurations: the generalization of the Jacobi sequence. Their detailed study will have to wait for other, 
possible numerical, studies. 



We have started the investigation of equilibrium figures of rotating fluids in higher dimensions. We used two 
appealingly simple models: one in which the fluid is held together by surface tension and the other in which self- 
gravity is the cohesive force. We have developed the basic tools to start an investigation of this important subject. 
Nevertheless, this work represents only a first step towards the understanding of the behavior of rotating fluids in 
higher dimensions. Many questions remain to be addressed. It is clear that the general behavior observed for D = 3 
still holds for higher D. It is natural to expect that (for instance for fluid drops), a non axisymmetric figure bifurcates 
from the point where the axisymmetric figure starts to be unstable, the point at which relation (|1C)4|) saturates. What 
is the shape of this figure? Is the instability point really a bifurcation point? This would allow one to generalize 
Figure |3 to higher dimensions. Furthermore, our study must be extended to encompass more general perturbations 
and more general angular momentum configurations. Another pressing issue is to determine whether there is always 
a toroidal figure of equilibrium in general D dimensions: could it be that for certain D a drop always encloses the 
origin regardless of its rotation? If so, this could be an excellent model for ultra-rotating Myers-Perry black holes, for 
which the angular momentum is unbounded. Concerning self-gravitating spheroids, there are many more unanswered 
questions. The formalism to study self-gravitating fluid drops must be further developed to find the actual equilibrium 
shapes as well as their stability. This would lead to a generalization of the MacLaurin and Jacobi sequences to higher 
dimensions. It would also be exciting to study numerically self-gravitating toroidal configurations (the analogue of 
the Dyson rings). 



The work presented here also allows one to better understand how gravity (even in the strong field regime) behaves 
in higher dimensions. Given the number of new solutions and instabilities discovered in this setting over the last 
few years, this is a particularly fascinating topic. A lesson to learn from the examples shown here (the Rayleigh- 
Plateau instability, the toroidal configurations, the DCF instability, the Dyson rings, etecetera) is that many attractive 
forces, like Newtonian gravity and surface tension, seem to mimic strong gravity situations very well. There is probably 
something more to it, when one shows that surface tension mimics so well the behavior of horizons, but this remains to 
be explained. A particularly promising proposal |65j . which gathers support from our work, is that the physics of fluids 
and black holes can be related semi-quantitativcly for smooth enough objects. If the object is very inhomogeneous 
or is highly distorted, then the gravitational self-interaction between different parts of the object will be strong, 
and then one expects that the details of the (fully non-linear) dynamics of the interaction will matter. So, for very 
inhomogeneous objects, general relativity or fluid dynamics would yield very different results. This means that the 
fluid approach won't, in principle, be useful for cases such as: the deeply non-linear evolution of the GL instability, 
when large inhomogeneities develop; the Myers-Perry black holes close to the extremal singularity (in the cases they 
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admit a Kerr bound); thick black rings, and in general black rings in the region where non- uniqueness occurs. This is 
unfortunate since these are very intriguing and important regimes. Nevertheless, a plethora of important and smooth 
objects (like the ones studied here) should then be described by these very simple models, and could even be a starting 
point for more complex situations. 

2. The bifurcation diagram of black objects 

In the Introduction, the suggestion is made that bifurcation diagram in Fig. ^ is m some sense universal, the 
qualitative behavior holding for most objects held together by some cohesive force. The natural generalization is to 
conjecture it may describe bifurcation diagrams of black objects. We would like to comment on the limitation of that 
generalization |79| . 

Kerr bound 

The first point to address concerns the existence of a Kerr bound. If the bifurcation diagram Fig. ^ were indeed 
universal, the angular momentum of black holes would be bounded. In fact, the results in Section Ivl imply that there 
is a bound on J in any of the cases I and II. We know however that the angular momentum of black holes in "case I" is 
unbounded (for D > 4) [H,^^ whereas in case II it is bounded [3^ (see for instance Eq. 7 in 68] ). Thus the surface 
tension model does not accurately describe ultra-spinning black holes, for the cases studied numerically (D = 3, 4, 5). 
Accurate numerical results are needed for higher dimensions, as well as results for self-gravitating fluids. This would 
clear the issue. 

Instabilities 

The second point concerns the instability of the solutions. The results available in the literature concerning black 
hole instability can be summarized as 

(i) For D > 4, case I, there should be both axisymmetric and non-axisymmetric instabilities 

(ii) For D = 4 spatial dimensions, case I, the analysis of ^3 does not apply; it was been suggested in ^(| that 
there could be an axisymmetric instability before the Kerr bound, ending in a black ring. 

(iii) For case II, the study in [bsT ] suggests no instability. 

The studies in [42| and focused only on a special kind of perturbation, related to a traceless tensor mode. This 
seems to belong to the same representation of our "toroidal mode" , even if one relates to the metric and the other to 
the fluid. 

In Section lYlII we found that a "case I" rotating drop is unstable. This agrees with the black hole picture (i) above. 
We were not able to extend the discussion to rotating drops in case II. So far both pictures (black hole and fluids) 
agree. The long-known discrepancy is the D = 3 case: Kerr black holes in D = 3 spatial dimensions are stable, 
whereas the fluid picture would predict instability. Another point of discrepancy is that the most general vacuum 
solution of Einstein's equations in 3 spatial dimensions is the Kerr family, with spherical topology so black rings are 
absent in this case. We can thus argue that D — 3 is a remarkable exception. 

Non-axisymmetric figures of equilibrium 

The final point we would like to touch upon is the question of existence of non-axisymmetric figures of equilibrium. 
In the simple fluid models (see Figs. ^ and [3J there is a family of non-axisymmetric stationary figures; one of the 
main points of our paper was to show that this holds for fluids in general higher dimensions, too. We went on to 
find non-axisymmetric instability which we argued brings the drop to a non-axisymmetric stationary configuration (in 
the case of gravitating fluid, we proved that there exist stationary, non-axisymmetric configurations analogous to the 
Jacobi sequence). However, it is shown in |69j that higher dimensional rotating black objects must be axisymmetric. 
Therefore stationary non-axisymmetric black solutions do not exist, and one would be led to think that the bifurcation 
diagram for fluids is completely misleading. 

The explanation for the disagreement is very simple, and lies in the fact that one is expecting from the model 
something it was not built to do: a more realistic analogy would have to take emission of gravitational waves into 
account, never considered up to now. Since non-axisymmetric objects radiate gravitational waves, we would not 
expect to find stationary equilibrium solutions of that kind. Indeed, for three spatial dimensions this problem was 
considered by Chandrasekhar |7dj| . Chandrasekhar proved that when one takes emission of gravitational waves into 
account (all other "curved spacetime" effects being neglected), the Jacobi configuration is unstable and evolves in 
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the direction of increasing angular velocity, approaching the bifurcation point. Furthermore, he proved that radiation 
reaction does not make the fluid unstable past the point of bifurcation. Thus the bifurcation diagram for black objects 
would be drawn by first deleting all the non-axisymmetric figures. 

It is reasonable to think that this result, i.e. the non-existence of non-axisymmetric stationary fluid configurations 
in general relativity, can be generalized to higher dimensions. If this is the case, this would be the analogue of the 
no-go result of 69] for higher dimensional rotating black holes. 



C. Future work 



The two examples we have studied (surface tension and self-gravitating fluids) can therefore improve and strengthen 
our understanding of black objects. The fact that rotating fluid drops develop instabilities at finite (usually small) 
rotation seems to suggest that black holes cannot have large rotation. This was also put forward by Emparan and 
Myers ^3] m connection with the Gregory-Laflamme instability. The examples shown here and the classical three- 
dimensional results discussed in the Introduction suggest that the final state of this instability is a non axi-symmetric 
figure. Some evidence was put forward already in |42| and this work adds some more evidence to that possibility. 
In fact, it might not even be possible to find matter configurations with high enough rotation to form these ultra- 
spinning black holes following collapse: we prove in Appendix [U] that self-gravitating, rotating homogeneous fluids 
always have a bounded angular velocity (this is a generalization of a classic result by Poincare). Future work along 
these lines (mimicking strong gravity situations with analogue models) could start by investigating the black string 
phase transition [52ll66ll6^ | within either the surface tension or the self-gravity models. Another potentially interesting 
topic would be to investigate the collapse of toroidal fluid configurations to form black rings. 
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APPENDIX A: OSCILLATION MODES OF A SPHERICAL DROP HELD TOGETHER BY SURFACE 

TENSION 



What are the oscillation modes of a spherical drop held together b y su rface tension? In three spatial dimensions 
the answer to this classical problem can be found for instance in |7ll l72| . We will now generalize it to an arbitrary 
number of dimensions. Take a spherical drop with radius Rq. Define again 

t = ^. (Ai) 

P 

where 5P is the change in pressure inside the drop (the outside pressure is taken to be zero), after it is perturbed by 
a displacement of the form 

r = R + ee wt r(angles) , (A2) 
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where K(angles) are hyper-spherical harmonics, with eigenvalue 1(1 + D — 2). The equations governing small pertur- 
bations are 

Po + SP = Tdivn, (A3) 

^ = -gradT, (A4) 

divu = 0, (A5) 
SP 

T = — . (A6) 
P 

With the displacement 1A2|) . equation (|A3|) can be solved to yield (for details we refer the reader to Landau and 
Lifshitz [72) near r — i?o and to first order in e 

SP T 

— = - lss -(!-i)(l + D- l)ee^Y (angles) . (A7) 

p RqP 



Therefore 

T = ' " ^Vee wt y (angles) . (A8) 



On the other hand, observing that V 2 =0we may immediately write 

5P 



eF p - r r(angles)e wt , (A9) 



with F p a constant to be determined. We thus have from the definition of T that 



Equation 1A4(I then yields 

1 d 



T = -^eF p y(angles)e wt . (A10) 



T , (All) 

u> ar 



Now, for this velocity to be compatible with the radial perturbation l|A2(l one must have 

Rou) 2 



F p = y- . (A12) 



Inserting this in (|A10|I and equating the resulting expression with (|A8|) we get 

u 2 = _W-W + D-D t (A13) 



which is a generalization of the three-dimensional result found in 0, |j3 ■ Note that the drop is always stable (lo 
is imaginary) and that I = and 1 modes do not exist. The first corresponds to spherically symmetric pulsation of 
the drop, which violates the constant volume assumption and the second corresponds to a translation of the drop as 
whole. 



APPENDIX B: OSCILLATION MODES OF A SELF-GRAVITATING, NON-ROTATING SPHERE 

Here we analyze the modes of a self-gravitating fluid sphere, in which case we can do without the virial formalism 
and find an exact analytical expression for the oscillation frequencies. We aim at generalizing the four-dimensional 
result by Kelvin and others [7l| . The analysis from Section IIXI carries over with minimal changes. 

The gravitational potential S(x) inside and outside the globe satisfy respectively 



V 2 Z?(x) 1 = -(D-2)C D G P , 
V 2 2?(x) e = 0. 



(Bl) 
(B2) 
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In the unperturbed state we get 



(D - 2)C D G P ( DRl 



«Wo = ^^- (B4) 



l(--BW) o. ,Br„ 



Hydrostatic equilibrium demands that the pressure Po satisfies 

dr \ p 

which yields, upon requiring Po(r = Rq) = 0, 

H- ffiz^ W-^). (B6, 

Suppose now the system is subjected to a perturbation, of general form 

r = R + ee^'y (angles) , (B7) 

where Y"(angles) are hyper-spherical harmonics, with eigenvalue 1(1 + D — 2). 
The equations governing small perturbations are 

(B8) 

IT , (B9) 

(BIO) 

5B(x) 1 . (Bll) 

Here u stands for the velocity of the perturbed fluid. Let's solve for the first of these equations. We find 

B(x) e = B(x)§ + eAr- l - D+2 Y (a,ngles)e ut , (B12) 
B(x)* = B(x)j + eSr'r(angles)e wt . (B13) 

The constants A and i? are determined by requiring continuity of the potential and its derivative on the boundary. 
By keeping only dominant terms in e we find 

_ (D-2)C D Gp x _, 
S - 2Z + £-2 ^ ' (B15) 



V 2 <5B(x) 


= o, 


9u 




~dt 


= -gi 


divu 


= o, 




SP 


T 






p 



Therefore we have 



W = e ^ r ^f^i?J- i r i r(angle S )e- . (B16) 



Observing that V ( 2£- ] = we may immediately write 



— = eF p - r F(angles)e wt , (B17) 
with F p a constant to be determined. We thus have from the definition of T, equation IjBllJ) . that 
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Equation i|B9|l then yields 

1 d 



T, (B19) 

oj dr 



Now, for this velocity to be compatible with the radial perturbation l|B7|) one must have 



Fp= (D-l)C D GpR _Ro^ 
p 21 + D -2 I y ' 



Imposing a vanishing pressure on the surface of the globe we have finally 



, 2{D - 2)C D Gp 1(1 - 1) 



"=-^5 -2!Td^2-> ^ 

which is a generalization of Kelvin's four-dimensional result |7l| to an arbitrary number of dimensions. As in the 
previous Section, the globe is always stable and I — and I = 1 modes do not exist. Result 1B21|I can be compared 
with the results for the modes of spherically symmetric black holes in higher dimensions |73i l74j . 



APPENDIX C: AN UPPER BOUND FOR THE ROTATION OF AN HOMOGENEOUS FLUID MASS 

In four dimensions there is an upper bound, derived by Poincare |4j for the maximum angular velocity of a rigidly 
rotating, homogeneous fluid mass in equilibrium. Consider a region V in three dimensional space bounded by a surface 
S. Let g^- denote differentiation of a function across S in the direction of the unit normal n e outward from V. The 

divergence theorem states that for any (sufficiently smooth) vector field W 

[ W.nd<r= I divWd 3 x. (CI) 
Js Jv 

If we choose W = gradS, where B is the gravitational potential, we get 

m ■ f 

v Jv 



,^ dcj= / V J Bd?x = - I 4irGpd 6 x, (C2) 



where we used Poisson's equation in the last equality. 

A necessary (but not sufficient) condition for equilibrium is that the force at every point in S be directed inward. 
But the total force is F r = J^- , where U is the modified potential 

U = B+Ul 2 ({x 1 ) 2 + (x 2 ) 2 ) (C3) 

(we choose a; 3 = z to be the axis of rotation). We then must require 

®L < => [ ^-da - / V 2 Ud 3 x = f (-AnGp + 2f! 2 ) d 3 x < . (C4) 
on e Js un e Jv Jv 

thus we have the Poincare's inequality [8(j 

O 2 < 2npG . (C5) 

This derivation can be easily generalized to D > 3 dimensional Euclidean space. If we assume the gravitational 
potential is (|175|l 



B(x)=Gjd D x' lx ^ D _ 2 , (C6) 

then Poisson's equation becomes 

V 2 B= -(D-2)C D Gp, (C7) 



3G 



where 



2tt d / 2 

Cd = , m , (C8) 



and Cor is the area of the (D — l)-sphere. Indeed, 



We now consider the two cases labeled by I and II, that is, 
I: one non-vanishing angular momentum; 
II: all angular momenta coincident, and D odd. 
The modified potential, in both cases, writes as 

u = B + ±n 2 Y / (x a ) 2 , (cio) 

a 

where a = 1, 2 in case I, a = 1, ...,£)— 1 in case II. Its Laplacian is 

V 2 W = + in case I, 

y 2 U = -C DP G + (D - 1)Q 2 in case II , (Cll) 
therefore Poincare's inequality becomes 

\l < ~^~P^ incase I, 

fl 2 < -jj^-jPG in case II . (C12) 

For angular velocities larger than this threshold, the gravitational pull is insufficient to balance the centrifugal force, 
and there body is disrupted. 

APPENDIX D: TOROIDAL PERTURBATIONS OF A ROTATING DROP 

Here we compute some useful quantities in the case of a rotating drop with one non-vanishing angular momentum, 
constant density p and with toroidal perturbations (following the denomination of Chandrasekhar |5|), i.e. 





= T ab x< 




= o, 


rj-iab 


/ a 







,)• (D1) 

with a, P constant parameters. We remind that, in the case with one non-vanishing angular momentum, the coordi- 
nates {x 1 } are decomposed as in (JTJ, l|BJ 

x l = (x a ,x A ); a = 1,2; A = 3,...,D (D2) 

and we have defined 

vu = ^(x^ + ix 2 ) 2 , 



z = (z 3 ) 2 + ■ ■ ■ + (x D ) 2 , (D3) 

so that 



x 1 — zd cos 9 , 



X 



2 = rosin 6. (D4) 
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The Z?-dimensional volume integration element is 

d D x — zuz D ~ 3 dmd9dzdQD-3 ■ (D5) 
If we integrate on the angular variables, we get 

d D x = 2nC D - 2 zuz D - 3 dwdz , (D6) 



S 1 xS 

where Cd-2 is the area of a unit (D — 3)-sphere l|C8l) . 

Let us now consider the (D — l)-dimensional surface S of equation z = f(zu). Since on S, dz — <pdw, the metric 
element on S is 

da 2 = (1 + <t>) 2 dw 2 + vj 2 d0 2 + [f(m)] 2 dn 2 D _ 3 . (D7) 
The surface integration element is then 

dS = a/1 + (f> 2 uj[f(zu)} D - 3 dmd0dn D ^ 3 . (D8) 

In particular, being 

Vi + 4> ™ 

we have 

dS a = n a dS = -(f>x a [f(m)} D - 3 dzud9dn D ^ 3 . (D10) 
After integration on the (D — 3)-sphere we get, 

/ dS a = -<f>x a [f(TXj)} D - 3 C D -2dmd6. (Dll) 
Js D - a 

Notice that in the D = 3 case we have dS a — —2<f>x a dmd0 where the factor 2 accounts for the two hemispheres of the 
drop. 



1. Virial integrals 

We can easily compute the virial integrals 

r ;J = / pCx J d D x. (D12) 
Jv 

From CDlJ it follows T AJ = 0; further more, 

/ x l x 3 d D x = if i (D13) 
Jv 

because the drop is symmetric with respect to reflections to any axis, therefore T a ' A = 0. 
We have 

/ (x 1 ) 2 d D x = / d9 cos 2 6 )C D -2 / dvjrn 3 / dzz°- 3 , 
Jv \Jq ) Jo Jo 

( '°~ 2 ^ dwvo 3 [f(w)] D - 2 = [ (x 2 ) 2 d D x , (D14) 



V 



D-2 

where w = a is the intersection of the drop with the x 1 — x 2 plane (namely, f(a) = 0); therefore 



r>l;l T^2;2 Cd-2 

I ' = —I ' = npa— 

H D-2 



/ w 3 [f{vj)] D ~ 2 dw, (D15) 
Jo 

r 1 ^ 2 = r 2 ^ 1 = £ vo 3 [f(m)] D ~ 2 dvo . (D16) 

We notice that these formulae are valid for every D > 3. When D > 3, Cd-2 represents the area of the unit 
(D — 3)-sphere; when D = 3, the factor C\ — 2 accounts for the two hemispheres composing the drop. 



2. Surface energy integrals 

Here we want to compute the quantities 5& 12 and 5& 11 — S& 22 where (see ffl\1 



S6 ab = _ T 



en\dS a + I x b S(n\ k )dS a - I x b n\ k ^-dS 



'dx a 



and 



dS a = -<j>x a [f(w)] u - 6 C D - 2 dmffl . 



Being £_ a = T ab x", we have that 



d0 



t-2 1 2 

iX - X dZ X 



2tt 



dS 



x — x — — r-X + X — — —X 



/ d6[T 2a x a x 1 — T la x 2 x a ] = J d6((x 1 ) 2 -(x 2 ) 2 )[3 = 0, 
Jo Jo 

■■ / d6[T la x a x 1 -T 2a x a x 2 - x l T cl x c + x 2 T c 
Jo 



dx 1 dx 2 

therefore the first and third integrals in (|D17|I cancel for the components 5& 12 and S& 11 — 56 22 



S& 12 = -T f 5{n k k )x 2 dS l , 
Js 

Se n -S6 22 = -T f S(n k k )x 1 dS 1 — T f 5(n k k )x 2 dS 2 
Js Js 



The perturbation of the divergence 8{n l 4 ) is given by (|51Jl 



= (5n%-^^ = (6n a ),a-T« b n% + (6n A U, 



++' : ! " -yr I / '" " " 



dx a 



[n n 1 \n — 1 n , 



5n A = [ n k n'^j ) n A = [n b n c T bc ]n A 



with 



We have 



then we find 



Therefore, 



where 



n a = -0- 



:^T+J 2 



tab 



(1 + 02)3/2 W 2 v ^02 W 3 ^/TT^VJ ' 



x a x b T ab 



0'(0 2 -2) 2 + 3(f) 2 D-3 



(1 + 2 ) 5 /2 + 02)3/2 z ( 1 + 02)3/2_ 

H n \i) = - [ a {( x1 ) 2 ~ i x2 f) + 2f3x 1 x 2 } Q(tu) = -w 2 (acos26 + f3sm26)Q(m), 

1 ^ 0'(0 2 -2) 2 + 302 £)-3 2 

(1+02)5/2 + ^ W (l + 02)3/2 "(1 + 02)3/5 
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The surface integrals 5& 12 , 5& 11 — S6 22 are then 

/•a pl-K 

6& 12 = TCd-2 dm d6(t)x 1 x 2 [f(m)} D - 3 5(n k k ), 
Jo Jo 

= -~PvC D - 2 f dm^Q^ifim)} - 3 , (D26) 
1 Jo 

5e 11 - 66 22 = TC D _ 2 I" dm / [(z 1 ) 2 - (x 2 ) 2 ] [f{m)} D - 3 S(n\ k ) , 

Jo Jo 

= -TmiC D -2 [ dm(j>m 4 Q(m)[f(m)] D - 3 . (D27) 
Jo 

APPENDIX E: NEWTON'S CONSTANT 

Here we would like to relate our adopted definition of Newton's constant to the one adopted in the by now classical 
work by Myers and Perry [3£|. They define a Newton's constant Gmp such that Einstein's equations read 

Riiv — -z9^vR — SttGhipT^ . (El) 

We will now derive "Newton's law" from this equation, in order to relate the different Newton constants. To derive 
Newton's law a = — VS(x), we consider the Newtonian limit of Einstein's equation (slowly moving particles, weak 
and static gravitational fields). 

A test particle follows the geodesies of the spacetimc 

d 2 r^ dx p d,r, a 

Under the "moving slowly" condition <C this is simplified to 



2 



d 2 x^ ,. ( dt\ 

r oo hr =0- (E3) 



dr 2 ' 00 V dr 
Due to the staticity assumption, 

Too = (9o9xo + dogox - d x g Q0 ) = - \g^d x g o , (E4) 
so that the geodesic equation simplifies to 

dr 2 2 3 "™ uu \drJ 
Lastly, the Newtonian limit implies weak gravitational fields, so we set 

9nv = Vfi"+f l nv, |/1|.k|<1. (E6) 
Thus we have (we can show that to first order g pv — rf v — h^) 



*" VW#Y = 0. (E5) 



(E7) 



and therefore 

dr 2 ~ 2'' ™ u \drj 
Dividing the spatial components of this equation by (^-) 2 we have finally 

d 2 x l 1 



dt 2 2 



d l h m , (E9) 
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which allows one to identify hoo = — 2£>(x). On the other hand, from Myers and Perry's equation (1.5), we get 

V 2 h oa = -167t^|g mf . (E10) 

Comparing our (|C7I) with (|E10|) . we get the relation between our adopted G and Gmp as (there's an extra minus sign 
that follows from the conventions adopted here for the relation force-potential) 

C D G = -j^-Gmp • (Ell) 
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